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ABSTRACT 

We set new constraints on a seven-dimensional space of cosmological parameters within the class of inflationary 
adiabatic models. We use the angular power spectrum of the cosmic microwave background measured over a wide 
range of £ in the first flight of the MAXIMA balloon-borne experiment (MAXIMA- 1) and the low t results from 
COBE/DMR. We find constraints on the total energy density of the universe, il = 1 .0+j] the physical density 
of baryons, il^h 2 = 0.03 ±0.01, the physical density of cold dark matter, fi c dm^ 2 = 0.2^ p and the spectral index 
of primordial scalar fluctuations, n s = 1.08 ±0.1, all at the 95% confidence level. By combining our results with 
measurements of high-redshift supernovae we constrain the value of the cosmological constant and the fractional 
amount of pressureless matter in the universe to 0.45 < 51a < 0.75 and 0.25 < il m < 0.50, at the 95% confidence 
level. Our results are consistent with a flat universe and the shape parameter deduced from large scale structure, 
and in marginal agreement with the baryon density from big bang nucleosynthesis. 

Subject headings: cosmic microwave background — cosmology: observations — large-scale structure of 
universe 
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1. INTRODUCTION of inflationary adiabatic models. 

™ , c • • , , This letter is structured as follows. In Section we de- 

The angular power spectrum of the cosmic microwave back- ., , . . r , . , a , 
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~LJ ground (CMB) temperature anisotropy depends strongly on , , 
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the values of most cosmological parameters and on the phys- , , • T o • ,, 
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Oical processes in the early universe; it is thus a powerful tool . . , ... J , J 

• , . ; At/ v ■ , , • „ MAXIMA- 1 and COBE/DMR data, to set constraints on the 

- — , to constrain cosmological models (see e.g., Kamionkowski & . „, , , . . . , , 
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O (see e.g., Lange et al. 2000, Tegmark & Zaldarriaga 2000, and ' , • , ^ 

e a. • \ sity in a cosmological constant, and pressureless matter, 
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» T tU . , „ , c . i , . il m . In Section W we summarize our results and discuss their 

rr, , In this letter we make use of the angular power spectrum . ... _ lj , , , , 

^4 a ■ iu c 4- a- u c a. \k a vth a u ii u implications for fundamental models of structure formation. 

■ measured in the first flight of the MAXIMA balloon-borne ex- F 

' periment, MAXIMA- 1 (described in a companion paper by 

• £h ! Hanany et al. 2000). This power spectrum covers a range 2. DATASET, METHODS AND PARAMETER SPACE 

>< 36 < £ < 785, which is the largest coverage in multipole space The MAX IMA-1 power spectrum is estimated in 10 bins 

H . from a single experiment to date. We include in our analysis spanning the range 36 < I < 785 (Hanany et al. 2000). In each 

... the 4-year COBE/DMR angular power spectrum (Gorski et al. bin ^ the spe ctrum is assigned a flat shape, £(£+l)Q/2ir = C B , 

1996), which spans the range 2 < £ < 34, to normalize the mod- whose amplitude is found by max i m i zing the likelihood of the 

els at large angular scales. We determine constraints on a seven- data using a quadratic estimator technique (Bond, Jaffe & Knox 

dimensional space of cosmological parameters within the class 199g) afj implemented in me MADCAP software package (Bor- 
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rill 1999). We find the most likely values of the cosmological 
parameters by maximizing the likelihood of the data. Follow- 
ing Bond, Jaffe & Knox (2000) we use a likelihood which is 
Gaussian in the quantity \n(C B +x B ), with x B related to the noise 
properties of the experiment. 

In this letter we only consider inflationary adiabatic models. 
There are two main reasons. First, the poor performance of 
alternative theories in reconciling measurements of clustering 
on galactic scales with measurements of CMB fluctuations on 
horizon scales has made inflation the favored contender as a 
fundamental theory of structure formation. Second, the angu- 
lar power spectra for inflationary models are easy to calculate 
for given cosmological parameters, especially since the advent 
of fast Einstein-Boltzmann solvers such as CMB FAST (Seljak 
& Zaldarriaga 1996) or CAMB (Lewis, Challinor & Lasenby 
2000). 

We consider a seven-dimensional space of parameters. This 
includes the amplitude of fluctuations at I = 10, Cm, the frac- 
tional densities of baryons, Sib, and cosmological constant, SIa, 
the total energy density of the universe, SI = Sl m + SI a (where 
the fractional density of pressureless matter is Sl m = Sib + ^cdm 
and Sl cdm is the fractional density of cold dark matter) the spec- 
tral index of primordial scalar fluctuations, n s , and the optical 
depth of reionization, t c . We use the following ranges and sam- 
pling: C m is continuous; SI = 0.1, 0.3, 0.5, 0.7, 0.8, 0.85, 0.9, 
0.95, 1, 1.05, 1.1, 1.15, 1.2, 1.3, 1.5; Sl b = 0.01, 0.03, 0.05, 
0.07, 0.09, 0.12, 0.20, 0.28; S1 A = 0.0, 0.1, 0.8, 0.9; n s = 
0.60, 0.64, • • •, 1.36, 1.40; r c = 0, 0.025, 0.05, 0.075, 0.1, 0.15, 
0.2, 0.3, 0.5. For the seventh parameter we have either used 
the Hubble parameter, H Q = 40, 45, • • •, 85, 90 km s _1 Mpc _1 , or 
the physical baryon density, Sl b /z 2 = 0.003125, 00625, 0.0125, 
0.0175, 0.02, 0.025, 0.030, 0.035, 0.04, 0.05, 0.075, 0.1, 0.15, 
0.2 (where h = Hq/(\QQ km s _1 Mpc -1 )). This approach allows 
us to check the importance of the parameters ranges, which are 
effectively used as a uniform prior in our analysis. When we 
use Hq as the free parameter we get relatively fine but restricted 
coverage of Hq and fl\,h 2 . When we use tt\,h 2 we consider a 
much more extended set of priors. To ensure that we are sam- 
pling the set of parameters with enough resolution, we include 
models corresponding to values of the parameters that are not 
on the grid by quadratically interpolating between power spec- 
tra. We have not allowed for the presence of massive neutrinos, 
since their effect on the angular power spectrum of the CMB 
would be negligible for the dataset used here (Dodelson, Gates 
& Stebbins 1996). 

We evaluate the likelihood for a subset of parameters by max- 
imizing it over all the remaining parameters. This approach, 
already adopted in previous analyses of this kind (Tegmark & 
Zaldarriaga 2000; Melchiorri et al. 1999), replaces the time- 
consuming integration that is required to marginalize over the 
unwanted parameters. In all our estimates we marginalize over 
the overall MAXIMA- 1 calibration uncertainty (Hanany et al. 
2000), assuming a Gaussian prior; we have also marginalized 
over the beam uncertainty reported in Hanany et al. (2000) but 
have found that it has a negligible effect on our results. 

3. CONSTRAINTS 

The position of the first peak in the angular power spec- 
trum of adiabatic models can be used to constrain the geometry 
of the universe (Doroshkevich, Zel'dovich & Sunyaev 1978). 
Features in the radiation pattern at recombination are set by 
the maximum distance sound waves have traveled at that time. 
The diameter-distance relation allows us to relate the physi- 



cal scale of these features to the angle they subtend on the 
sky and strongly depends on the curvature of the universe, 
which is determined by the total energy density of the universe 
SI = Sl m + SlA. Other cosmological parameters have much less 
effect on this relation. 




Fig. 1. — Likelihood function of the total energy density of the uni- 
verse Q. The solid line was obtained by maximizing over all other 
parameters over the ranges described in the text, while the dashed 
line was obtained by constraining Q,bh 2 = 0.0190 ± 0.0024 and Ho = 
65 ± 7 km s _1 Mpc~' . The intersections with the horizontal line give 
the bounds for 95% confidence. 



Figure 1 shows the likelihood of the total energy density 
of the universe, SI. We obtain the solid line if we maxi- 
mize over all the remaining parameters; we get the dashed line 
if we impose the big bang nucleosynthesis (BBN) constraint 
n b h z = 0.0190 ± 0.0024 (Buries et al. 1999; see also, Tytler et 
al. 2000), and restrict H = 65 ± 7 km s" 1 Mpc" 1 (Freedman 
1999). From the solid line we see that 0.7 < SI < 1.15 at the 
95% confidence level. 

We identified the other tightly constrained parameters from 
a principal component analysis of the likelihood function (Efs- 
tathiou & Bond 1999). In Figure 2 we plot the likelihoods for 
three well-constrained parameters, the physical baryon density, 
Slb^ 2 , the physical cold dark matter density, Sl c( j m /z 2 , and the 
spectral index of primordial scalar fluctuations, n s . The likeli- 
hood for each parameter was obtained by maximizing over all 
the remaining parameters. 

We find 0.02 < fl h h 2 < 0.04, 0.1 < Sl cdm /i 2 < 0.4, and 0.98 < 
n s < 1.18, all at the 95% confidence level. The constraint on 
Slb^ 2 is independent of the one coming from BBN (Buries et 
al. 1999) and agrees with it at about the 2a level. However, the 
most likely value of Sl b /z 2 emerging from our dataset is notice- 
ably higher than the one from BBN. Imposing the BBN prior to 
our dataset has the effect of moving the best fit region towards 
open models (see Figure 1). The constraint on n s indicates that 
the combined MAXIMA- 1 and COBE/DMR data are consistent 
with a scale invariant spectrum, with a smaller variance than 
that estimated using previous datasets (see e.g., Lange et al. 
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2000, Tegmark & Zaldarriaga 2000). This result should how- 
ever be interpreted with some caution, because it is sensitive 
to assumptions on the treatment of calibration uncertainties and 
depends in part on the absence of tensor modes in our analysis. 
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FIG. 2. — Likelihoods of three cosmological parameters: from top 
to bottom, the physical baryon density, Q.bh 2 , the physical cold dark 
matter density, Q c <i m h , and the scalar spectral index, n s . For each of 
these parameters, the likelihood was maximized over all the remaining 
parameters. The vertical band in the top panel represents the BBN con- 
straint fUi 2 = 0.0190 ± 0.0024 (Buries et al. 1999). The intersections 
with the horizontal line give the bounds for 95% confidence. 



One of the intriguing developments arising in contemporary 
observational cosmology is the possible existence of a smooth, 
negative-pressure component, for example a cosmological con- 
stant, driving an accelerated expansion of the universe (Perl- 
mutter et al. 1999; Riess et al. 1998). While CMB data are 
quite powerful in constraining the total energy density of the 
universe, a precise independent determination of tt m and 
is limited by the presence of the geometrical degeneracy dis- 
cussed in Efstathiou & Bond (1999) (which can be broken if a 
gravitational lensing imprint on the CMB is detected, see e.g., 
Stompor & Efstathiou 1999). We can however find the locus 
of models in the fi m — fl\ plane which are favored by the com- 
bined MAXIMA- 1 and COBE/DMR data. We do so in Fig- 
ure 3, where once again we maximize over all the remaining pa- 
rameters. The fact that the MAXIMA- 1 dataset extends to large 
values of £ helps to narrow the contours along the degeneracy 
direction 51 m + J7a = constant. Although a large portion of the 
fl m — £7a plane is ruled out, one clearly needs additional data to 
strongly break the degeneracy. We then overlay our results on 
the likelihood contours from the high-redshift supernovae data 
of Perlmutter et al. (1999) and Riess et al. (1998) and calculate 



the combined likelihood. The resulting constraints are: 0.45 < 
fi A < 0.75 and 0.25 < ft m < 0.50, at the 95% confidence level. 
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Fig. 3. — Constraints in the f2 m -fiA plane from the combined 
MAXIMA- 1 and COBE/DMR datasets. The likelihood at each point 
in the plane is taken by maximizing over the remaining five parame- 
ters. The borders of the shaded regions correspond to 0.32, 0.05 and 
0.01 of the peak value of the likelihood. These are overlaid on the 
bounds obtained from high redshift supernovae data (Perlmutter et al. 
1999; Riess et al. 1998). The closed contours are the confidence levels 
from the combined likelihood. The diagonal line corresponds to flat 
models. 



4. CONCLUSIONS 

Using the broad £-space coverage of the MAXIMA- 1 and 
COBE/DMR data we set constraints on a number of cosmologi- 
cal parameters within the class of inflationary adiabatic models. 
We find: n = 1.0^, n h h 2 = 0.03 ±0.01, n cim h 2 = 0.2^ and 
n s = 1.08 ±0.1 at the 95% confidence level. The constraints 
we have obtained on 57 are consistent with those of other re- 
cent analyses (Melchiorri et al. 1999; Dodelson & Knox 1999; 
Tegmark & Zaldarriaga 2000; Lange et al. 2000) and pro- 
vide further support to a flat universe. Combining our con- 
straints with those coming from measurements of supernovae at 
high redshift strongly favors a non-zero cosmological constant, 
17a = 0.60 ±0.15, and a pressureless matter density parameter 
fl m = 0.375 ± 0.125, both at the 95% confidence level. 

Our analysis was restricted to inflationary adiabatic mod- 
els. Several hints suggest that this may indeed be the right 
paradigm. The presence of a localized peak in our data in the 
region 150 < £ < 250 is in itself an evidence in favor of in- 
flationary adiabatic models; alternative theories either predict 
a broader peak at higher £ or a broad shelf at £ < 200 (see 
e.g., Knox & Page 1999). Furthermore, from a goodness-of-fit 
analysis (with the caveat that our data are not perfectly well- 
described by a Gaussian likelihood) we find that the class of 
models we have considered are perfectly consistent with the 
data. Our best fit has x 2= 41 using the 38 data points of the 
combined MAXIMA-1/COBE dataset (x 2 =7 using only the 10 
data points from MAXIMA-1). Within inflationary models, the 
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Standard CDM model and the Open CDM model with = 0.3 
both provide a poor fit to our data, having, respectively, x 2 =67 
(25) and x 2 =224 (159). The A-CDM 'concordance model' (Os- 
triker & Steinhardt 1995; Krauss & Turner 1995) has x 2 =50 (9). 
These results are summarized in Figure 4, where we plot the 
MAXIMA- 1 data overlaid with various cosmological models. 
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Fig. 4. — Angular power spectrum from MAXIMA- 1 (Hanany et 
al. 2000). The solid line is the best fit model of our analysis, having 
(A,ficdm,ftA, n s ,h) = (0.1,0.6,0.3, 1.08,0.53). The long-dashed line 
is the Standard CDM model (0.05,0.95,0,1,0.5). The short-dashed line 
is the A-CDM model with (0.05,0.35,0.6,1,0.65). The dot-dashed line 
is the Open CDM model with ft = 0.3. 



Finally, the agreement of our results with those obtained from 



other observations lends further credence to the possibility that 
we are considering the correct family of cosmological mod- 
els. First, we find an independent constraint on the physical 
baryon density of the universe: 0.02 < n b h 2 < 0.04, which, al- 
though favoring higher values than the BBN estimate of Buries 
et al. (1999), agrees with it within 2a. Second, we are able 
to constrain the amount of cold dark matter in the universe: 
^cdm^ 2 = 0.2^q 2 . If we combine this value with our constraint 
on n h h 2 and the MAXIMA- 1/COBE/SN la best fit value of 
fi m ~ 0.375 derived from Figure 3, we find that the "shape" 
parameter (commonly used in analysis of large scale structure; 
see e.g., Sugiyama 1995) is T = f2 m /z exp (-fib - fib/fim) = 
0.24 ±0.1. Together with our best estimate of ris = 1.08 ±0.1, 
this is consistent with the value T = 0.23-0.28(1 - which 
emerges from a completely independent analysis of galaxy cat- 
alogs (Viana & Liddle 1999; Liddle et al. 1995; Peacock & 
Dodds 1994). 
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ERRATUM 



In the Letter "Constraints on Cosmological Parameters from MAXIMA-1" by A. Balbi et al. (ApJ 545, LI [2001]), an error was 
found in the results: while it is claimed in the Letter that the optical depth of reionization, t c , took values bewteen and 0.5, an error 
in the computer code led to the results being based on t c = 0. We have regenerated the database of models with the correct (published) 
values of r c , namely, t c = 0, 0.025, 0.05, 0.075, 0.1, 0.15, 0.2, 0.3, 0.5. The best-fit model quoted in the Letter remains unchanged. 
However, we have found that the 95% confidence level constraints on the total energy density of the universe, the physical density 
of baryons and the scalar spectral index of primordial fluctuations change to: O = l-O^Jo, ^bh 2 = 0.030^[J;o[o and n s = 1.08^q o|, 
respectively. The likelihoods for the remaining parameters are effectively unchanged. In Figures 1 and 2, we present the figures 
which should replace Figures 1 and 2 of the original Letter. 
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Fig. 1. — Likelihood function of the total energy density of the uni- 
verse Q. The solid line was obtained by maximizing over all other 
parameters over the ranges described in the text, while the dashed 
line was obtained by constraining Q b h 2 = 0.0190 ± 0.0024 and H = 
65 ± 7 km s -1 Mpc~' . The intersections with the horizontal line give 
the bounds for 95% confidence. 
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FIG. 2. — Likelihoods of three cosmological parameters: from top 
to bottom, the physical baryon density, Q. b h 2 , the physical cold dark 
matter density, Q c a m h 2 , and the scalar spectral index, n s . For each of 
these parameters, the likelihood was maximized over all the remaining 
parameters. The vertical band in the top panel represents the BBN con- 
straint Q. h h 2 = 0.0190 ±0.0024 (Buries et al. 1999). The intersections 
with the horizontal line give the bounds for 95% confidence. 



